aenmd: annotating escape from nonsense-mediated decay for transcripts with protein-truncating variants.

SUMMARY: DNA changes that cause premature termination codons (PTCs) represent a large fraction of clinically relevant pathogenic genomic variation. Typically, PTCs induce transcript degradation by nonsense-mediated mRNA decay (NMD) and render such changes loss-of-function alleles. However, certain PTC-containing transcripts escape NMD and can exert dominant-negative or gain-of-function (DN/GOF) effects. Therefore, systematic identification of human PTC-causing variants and their susceptibility to NMD contributes to the investigation of the role of DN/GOF alleles in human disease. Here we present aenmd, a software for annotating PTC-containing transcript-variant pairs for predicted escape from NMD. aenmd is user-friendly and self-contained. It offers functionality not currently available in other methods and is based on established and experimentally validated rules for NMD escape; the software is designed to work at scale, and to integrate seamlessly with existing analysis workflows. We applied aenmd to variants in the gnomAD, Clinvar, and GWAS catalog databases and report the prevalence of human PTC-causing variants in these databases, and the subset of these variants that could exert DN/GOF effects via NMD escape. AVAILABILITY AND IMPLEMENTATION: aenmd is implemented in the R programming language. Code is available on GitHub as an R-package (github.com/kostkalab/aenmd.git), and as a containerized command-line interface (github.com/kostkalab/aenmd_cli.git).

An intrinsically interpretable neural network architecture for sequence-to-function learning.

MOTIVATION: Sequence-based deep learning approaches have been shown to predict a multitude of functional genomic readouts, including regions of open chromatin and RNA expression of genes. However, a major limitation of current methods is that model interpretation relies on computationally demanding post hoc analyses, and even then, one can often not explain the internal mechanics of highly parameterized models. Here, we introduce a deep learning architecture called totally interpretable sequence-to-function model (tiSFM). tiSFM improves upon the performance of standard multilayer convolutional models while using fewer parameters. Additionally, while tiSFM is itself technically a multilayer neural network, internal model parameters are intrinsically interpretable in terms of relevant sequence motifs. RESULTS: We analyze published open chromatin measurements across hematopoietic lineage cell-types and demonstrate that tiSFM outperforms a state-of-the-art convolutional neural network model custom-tailored to this dataset. We also show that it correctly identifies context-specific activities of transcription factors with known roles in hematopoietic differentiation, including Pax5 and Ebf1 for B-cells, and Rorc for innate lymphoid cells. tiSFM's model parameters have biologically meaningful interpretations, and we show the utility of our approach on a complex task of predicting the change in epigenetic state as a function of developmental transition. AVAILABILITY AND IMPLEMENTATION: The source code, including scripts for the analysis of key findings, can be found at https://github.com/boooooogey/ATAConv, implemented in Python.

Spatial transcriptome profiling uncovers metabolic regulation of left-right patterning.

Left-right patterning disturbance can cause severe birth defects, but it remains least understood of the three body axes. We uncovered an unexpected role for metabolic regulation in left-right patterning. Analysis of the first spatial transcriptome profile of left-right patterning revealed global activation of glycolysis, accompanied by right-sided expression of Bmp7 and genes regulating insulin growth factor signaling. Cardiomyocyte differentiation was left-biased, which may underlie the specification of heart looping orientation. This is consistent with known Bmp7 stimulation of glycolysis and glycolysis suppression of cardiomyocyte differentiation. Liver/lung laterality may be specified via similar metabolic regulation of endoderm differentiation. Myo1d , found to be left-sided, was shown to regulate gut looping in mice, zebrafish, and human. Together these findings indicate metabolic regulation of left-right patterning. This could underlie high incidence of heterotaxy-related birth defects in maternal diabetes, and the association of PFKP, allosteric enzyme regulating glycolysis, with heterotaxy. This transcriptome dataset will be invaluable for interrogating birth defects involving laterality disturbance.

aenmd: Annotating escape from nonsense-mediated decay for transcripts with protein-truncating variants.

DNA changes that cause premature termination codons (PTCs) represent a large fraction of clinically relevant pathogenic genomic variation. Typically, PTCs induce a transcript's degradation by nonsense-mediated mRNA decay (NMD) and render such changes loss-of-function alleles. However, certain PTC-containing transcripts escape NMD and can exert dominant-negative or gain-of-function (DN/GOF) effects. Therefore, systematic identification of human PTC-causing variants and their susceptibility to NMD contributes to the investigation of the role of DN/GOF alleles in human disease. Here we present aenmd, a software for annotating PTC-containing transcript-variant pairs for predicted escape from NMD. aenmd is user-friendly and self-contained. It offers functionality not currently available in other methods and is based on established and experimentally validated rules for NMD escape; the software is designed to work at scale, and to integrate seamlessly with existing analysis workflows. We applied aenmd to variants in the gnomAD, Clinvar, and GWAS catalog databases and report the prevalence of human PTC-causing variants in these databases, and the subset of these that could exert DN/GOF effects via NMD escape. Availability and implementation: aenmd is implemented in the R programming language. Code is available on GitHub as an R package (github.com/kostkalab/aenmd.git), and as a containerized command-line interface (github.com/kostkalab/aenmd_cli.git).

An intrinsically interpretable neural network architecture for sequence to function learning.

MOTIVATION: Sequence-based deep learning approaches have been shown to predict a multitude of functional genomic readouts, including regions of open chromatin and RNA expression of genes. However, a major limitation of current methods is that model interpretation relies on computationally demanding post hoc analyses, and even then, one can often not explain the internal mechanics of highly parameterized models. Here, we introduce a deep learning architecture called tiSFM (totally interpretable sequence to function model). tiSFM improves upon the performance of standard multi-layer convolutional models while using fewer parameters. Additionally, while tiSFM is itself technically a multi-layer neural network, internal model parameters are intrinsically interpretable in terms of relevant sequence motifs. RESULTS: We analyze published open chromatin measurements across hematopoietic lineage cell-types and demonstrate that tiSFM outperforms a state-of-the-art convolutional neural network model custom-tailored to this dataset. We also show that it correctly identifies context specific activities of transcription factors with known roles in hematopoietic differentiation, including Pax5 and Ebf1 for B-cells, and Rorc for innate lymphoid cells. tiSFM's model parameters have biologically meaningful interpretations, and we show the utility of our approach on a complex task of predicting the change in epigenetic state as a function of developmental transition. AVAILABILITY AND IMPLEMENTATION: The source code, including scripts for the analysis of key findings, can be found at https://github.com/boooooogey/ATAConv, implemented in Python.

Vaeda computationally annotates doublets in single-cell RNA sequencing data.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) continues to expand our knowledge by facilitating the study of transcriptional heterogeneity at the level of single cells. Despite this technology's utility and success in biomedical research, technical artifacts are present in scRNA-seq data. Doublets/multiplets are a type of artifact that occurs when two or more cells are tagged by the same barcode, and therefore they appear as a single cell. Because this introduces non-existent transcriptional profiles, doublets can bias and mislead downstream analysis. To address this limitation, computational methods to annotate and remove doublets form scRNA-seq datasets are needed. RESULTS: We introduce vaeda (Variational Auto-Encoder for Doublet Annotation), a new approach for computational annotation of doublets in scRNA-seq data. Vaeda integrates a variational auto-encoder and Positive-Unlabeled learning to produce doublet scores and binary doublet calls. We apply vaeda, along with seven existing doublet annotation methods, to 16 benchmark datasets and find that vaeda performs competitively in terms of doublet scores and doublet calls. Notably, vaeda outperforms other python-based methods for doublet annotation. Altogether, vaeda is a robust and competitive method for scRNA-seq doublet annotation and may be of particular interest in the context of python-based workflows. AVAILABILITY AND IMPLEMENTATION: Vaeda is available at https://github.com/kostkalab/vaeda, and the version used for the results we present here is archived at zenodo (https://doi.org/10.5281/zenodo.7199783). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Single-cell transcriptomic analysis identifies murine heart molecular features at embryonic and neonatal stages.

Heart development is a continuous process involving significant remodeling during embryogenesis and neonatal stages. To date, several groups have used single-cell sequencing to characterize the heart transcriptomes but failed to capture the progression of heart development at most stages. This has left gaps in understanding the contribution of each cell type across cardiac development. Here, we report the transcriptional profile of the murine heart from early embryogenesis to late neonatal stages. Through further analysis of this dataset, we identify several transcriptional features. We identify gene expression modules enriched at early embryonic and neonatal stages; multiple cell types in the left and right atriums are transcriptionally distinct at neonatal stages; many congenital heart defect-associated genes have cell type-specific expression; stage-unique ligand-receptor interactions are mostly between epicardial cells and other cell types at neonatal stages; and mutants of epicardium-expressed genes Wt1 and Tbx18 have different heart defects. Assessment of this dataset serves as an invaluable source of information for studies of heart development.

Non-negative Independent Factor Analysis disentangles discrete and continuous sources of variation in scRNA-seq data.

MOTIVATION: Single-cell RNA-seq analysis has emerged as a powerful tool for understanding inter-cellular heterogeneity. Due to the inherent noise of the data, computational techniques often rely on dimensionality reduction (DR) as both a pre-processing step and an analysis tool. Ideally, DR should preserve the biological information while discarding the noise. However, if the DR is to be used directly to gain biological insight it must also be interpretable-that is the individual dimensions of the reduction should correspond to specific biological variables such as cell-type identity or pathway activity. Maximizing biological interpretability necessitates making assumption about the data structures and the choice of the model is critical. RESULTS: We present a new probabilistic single-cell factor analysis model, Non-negative Independent Factor Analysis (NIFA), that incorporates different interpretability inducing assumptions into a single modeling framework. The key advantage of our NIFA model is that it simultaneously models uni- and multi-modal latent factors, and thus isolates discrete cell-type identity and continuous pathway activity into separate components. We apply our approach to a range of datasets where cell-type identity is known, and we show that NIFA-derived factors outperform results from ICA, PCA, NMF and scCoGAPS (an NMF method designed for single-cell data) in terms of disentangling biological sources of variation. Studying an immunotherapy dataset in detail, we show that NIFA is able to reproduce and refine previous findings in a single analysis framework and enables the discovery of new clinically relevant cell states. AVAILABILITY AND IMPLEMENTATION: NFIA is a R package which is freely available at GitHub (https://github.com/wgmao/NIFA). The test dataset is archived at https://zenodo.org/record/6286646. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational profiling of hiPSC-derived heart organoids reveals chamber defects associated with NKX2-5 deficiency.

Heart organoids have the potential to generate primary heart-like anatomical structures and hold great promise as in vitro models for cardiac disease. However, their properties have not yet been fully studied, which hinders their wide spread application. Here we report the development of differentiation systems for ventricular and atrial heart organoids, enabling the study of heart diseases with chamber defects. We show that our systems generate chamber-specific organoids comprising of the major cardiac cell types, and we use single cell RNA sequencing together with sample multiplexing to characterize the cells we generate. To that end, we developed a machine learning label transfer approach leveraging cell type, chamber, and laterality annotations available for primary human fetal heart cells. We then used this model to analyze organoid cells from an isogeneic line carrying an Ebstein's anomaly associated genetic variant in NKX2-5, and we successfully recapitulated the disease's atrialized ventricular defects. In summary, we have established a workflow integrating heart organoids and computational analysis to model heart development in normal and disease states.

Uncompensated mitochondrial oxidative stress underlies heart failure in an iPSC-derived model of congenital heart disease.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease with 30% mortality from heart failure (HF) in the first year of life, but the cause of early HF remains unknown. Induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CM) from patients with HLHS showed that early HF is associated with increased apoptosis, mitochondrial respiration defects, and redox stress from abnormal mitochondrial permeability transition pore (mPTP) opening and failed antioxidant response. In contrast, iPSC-CM from patients without early HF showed normal respiration with elevated antioxidant response. Single-cell transcriptomics confirmed that early HF is associated with mitochondrial dysfunction accompanied with endoplasmic reticulum (ER) stress. These findings indicate that uncompensated oxidative stress underlies early HF in HLHS. Importantly, mitochondrial respiration defects, oxidative stress, and apoptosis were rescued by treatment with sildenafil to inhibit mPTP opening or TUDCA to suppress ER stress. Together these findings point to the potential use of patient iPSC-CM for modeling clinical heart failure and the development of therapeutics.

Chromatin accessibility and microRNA expression in nephron progenitor cells during kidney development.

Mammalian nephrons originate from a population of nephron progenitor cells, and changes in these cells' transcriptomes contribute to the cessation of nephrogenesis, an important determinant of nephron number. To characterize microRNA (miRNA) expression and identify putative cis-regulatory regions, we collected nephron progenitor cells from mouse kidneys at embryonic day 14.5 and postnatal day zero and assayed small RNA expression and transposase-accessible chromatin. We detect expression of 1104 miRNA (114 with expression changes), and 46,374 chromatin accessible regions (2103 with changes in accessibility). Genome-wide, our data highlight processes like cellular differentiation, cell migration, extracellular matrix interactions, and developmental signaling pathways. Furthermore, they identify new candidate cis-regulatory elements for Eya1 and Pax8, both genes with a role in nephron progenitor cell differentiation. Finally, we associate expression-changing miRNAs, including let-7-5p, miR-125b-5p, miR-181a-2-3p, and miR-9-3p, with candidate cis-regulatory elements and target genes. These analyses highlight new putative cis-regulatory loci for miRNA in nephron progenitors.

Single-cell RNA sequencing reveals differential cell cycle activity in key cell populations during nephrogenesis.

The kidney is a complex organ composed of more than 30 terminally differentiated cell types that all are required to perform its numerous homeostatic functions. Defects in kidney development are a significant cause of chronic kidney disease in children, which can lead to kidney failure that can only be treated by transplant or dialysis. A better understanding of molecular mechanisms that drive kidney development is important for designing strategies to enhance renal repair and regeneration. In this study, we profiled gene expression in the developing mouse kidney at embryonic day 14.5 at single-cell resolution. Consistent with previous studies, clusters with distinct transcriptional signatures clearly identify major compartments and cell types of the developing kidney. Cell cycle activity distinguishes between the "primed" and "self-renewing" sub-populations of nephron progenitors, with increased expression of the cell cycle-related genes Birc5, Cdca3, Smc2 and Smc4 in "primed" nephron progenitors. In addition, augmented expression of cell cycle related genes Birc5, Cks2, Ccnb1, Ccnd1 and Tuba1a/b was detected in immature distal tubules, suggesting cell cycle regulation may be required for early events of nephron patterning and tubular fusion between the distal nephron and collecting duct epithelia.

Apoptosis of hematopoietic progenitor-derived adipose tissue-resident macrophages contributes to insulin resistance after myocardial infarction.

Patients with insulin resistance have high risk of cardiovascular disease such as myocardial infarction (MI). However, it is not known whether MI can initiate or aggravate insulin resistance. We observed that patients with ST-elevation MI and mice with MI had de novo hyperglycemia and features of insulin resistance, respectively. In mouse models of both myocardial and skeletal muscle injury, we observed that the number of visceral adipose tissue (VAT)-resident macrophages decreased because of apoptosis after these distant organ injuries. Patients displayed a similar decrease in VAT-resident macrophage numbers and developed systemic insulin resistance after ST-elevation MI. Loss of VAT-resident macrophages after MI injury led to systemic insulin resistance in non-diabetic mice. Danger signaling-associated protein high mobility group box 1 was released by the dead myocardium after MI in rodents and triggered macrophage apoptosis via Toll-like receptor 4. The VAT-resident macrophage population in the steady state in mice was transcriptomically distinct from macrophages in the brain, skin, kidney, bone marrow, lungs, and liver and was derived from hematopoietic progenitor cells just after birth. Mechanistically, VAT-resident macrophage apoptosis and de novo insulin resistance in mouse models of MI were linked to diminished concentrations of macrophage colony-stimulating factor and adiponectin. Collectively, these findings demonstrate a previously unappreciated role of adipose tissue-resident macrophages in sensing remote organ injury and promoting MI pathogenesis.

Deletion of hypoxia-responsive microRNA-210 results in a sex-specific decrease in nephron number.

Low nephron number results in an increased risk of developing hypertension and chronic kidney disease. Intrauterine growth restriction is associated with a nephron deficit in humans, and is commonly caused by placental insufficiency, which results in fetal hypoxia. The underlying mechanisms by which hypoxia impacts kidney development are poorly understood. microRNA-210 is the most consistently induced microRNA in hypoxia and is known to promote cell survival in a hypoxic environment. In this study, the role of microRNA-210 in kidney development was evaluated using a global microRNA-210 knockout mouse. A male-specific 35% nephron deficit in microRNA-210 knockout mice was observed. Wnt/ß-catenin signaling, a pathway crucial for nephron differentiation, was misregulated in male kidneys with increased expression of the canonical Wnt target lymphoid enhancer binding factor 1. This coincided with increased expression of caspase-8-associated protein 2, a known microRNA-210 target and apoptosis signal transducer. Together, these data are consistent with a sex-specific requirement for microRNA-210 in kidney development.

scds: computational annotation of doublets in single-cell RNA sequencing data.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) technologies enable the study of transcriptional heterogeneity at the resolution of individual cells and have an increasing impact on biomedical research. However, it is known that these methods sometimes wrongly consider two or more cells as single cells, and that a number of so-called doublets is present in the output of such experiments. Treating doublets as single cells in downstream analyses can severely bias a study's conclusions, and therefore computational strategies for the identification of doublets are needed. RESULTS: With scds, we propose two new approaches for in silico doublet identification: Co-expression based doublet scoring (cxds) and binary classification based doublet scoring (bcds). The co-expression based approach, cxds, utilizes binarized (absence/presence) gene expression data and, employing a binomial model for the co-expression of pairs of genes, yields interpretable doublet annotations. bcds, on the other hand, uses a binary classification approach to discriminate artificial doublets from original data. We apply our methods and existing computational doublet identification approaches to four datasets with experimental doublet annotations and find that our methods perform at least as well as the state of the art, at comparably little computational cost. We observe appreciable differences between methods and across datasets and that no approach dominates all others. In summary, scds presents a scalable, competitive approach that allows for doublet annotation of datasets with thousands of cells in a matter of seconds. AVAILABILITY AND IMPLEMENTATION: scds is implemented as a Bioconductor R package (doi: 10.18129/B9.bioc.scds). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

FSTL-1 Attenuation Causes Spontaneous Smoke-Resistant Pulmonary Emphysema.

Rationale: The role of FSTL-1 (follistatin-like 1) in lung homeostasis is unknown.Objectives: We aimed to define the impact of FSTL-1 attenuation on lung structure and function and to identify FSTL-1-regulated transcriptional pathways in the lung. Further, we aimed to analyze the association of FSTL-1 SNPs with lung disease.Methods: FSTL-1 hypomorphic (FSTL-1 Hypo) mice underwent lung morphometry, pulmonary function testing, and micro-computed tomography. Fstl1 expression was determined in wild-type lung cell populations from three independent research groups. RNA sequencing of wild-type and FSTL-1 Hypo mice identified FSTL-1-regulated gene expression, followed by validation and mechanistic in vitro examination. FSTL1 SNP analysis was performed in the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) cohort.Measurements and Main Results: FSTL-1 Hypo mice developed spontaneous emphysema, independent of smoke exposure. Fstl1 is highly expressed in the lung by mesenchymal and endothelial cells but not immune cells. RNA sequencing of whole lung identified 33 FSTL-1-regulated genes, including Nr4a1, an orphan nuclear hormone receptor that negatively regulates NF-κB (nuclear factor-κB) signaling. In vitro, recombinant FSTL-1 treatment of macrophages attenuated NF-κB p65 phosphorylation in an Nr4a1-dependent manner. Within the COPDGene cohort, several SNPs in the FSTL1 region corresponded to chronic obstructive pulmonary disease and lung function.Conclusions: This work identifies a novel role for FSTL-1 protecting against emphysema development independent of smoke exposure. This FSTL-1-deficient emphysema implicates regulation of immune tolerance in lung macrophages through Nr4a1. Further study of the mechanisms involving FSTL-1 in lung homeostasis, immune regulation, and NF-κB signaling may provide additional insight into the pathophysiology of emphysema and inflammatory lung diseases.

Control of cytokinesis by ß-adrenergic receptors indicates an approach for regulating cardiomyocyte endowment.

One million patients with congenital heart disease (CHD) live in the United States. They have a lifelong risk of developing heart failure. Current concepts do not sufficiently address mechanisms of heart failure development specifically for these patients. Here, analysis of heart tissue from an infant with tetralogy of Fallot with pulmonary stenosis (ToF/PS) labeled with isotope-tagged thymidine demonstrated that cardiomyocyte cytokinesis failure is increased in this common form of CHD. We used single-cell transcriptional profiling to discover that the underlying mechanism of cytokinesis failure is repression of the cytokinesis gene ECT2, downstream of ß-adrenergic receptors (ß-ARs). Inactivation of the ß-AR genes and administration of the ß-blocker propranolol increased cardiomyocyte division in neonatal mice, which increased the number of cardiomyocytes (endowment) and conferred benefit after myocardial infarction in adults. Propranolol enabled the division of ToF/PS cardiomyocytes in vitro. These results suggest that ß-blockers could be evaluated for increasing cardiomyocyte division in patients with ToF/PS and other types of CHD.

Large-scale inference of competing endogenous RNA networks with sparse partial correlation.

MOTIVATION: MicroRNAs (miRNAs) are important non-coding post-transcriptional regulators that are involved in many biological processes and human diseases. Individual miRNAs may regulate hundreds of genes, giving rise to a complex gene regulatory network in which transcripts carrying miRNA binding sites act as competing endogenous RNAs (ceRNAs). Several methods for the analysis of ceRNA interactions exist, but these do often not adjust for statistical confounders or address the problem that more than one miRNA interacts with a target transcript. RESULTS: We present SPONGE, a method for the fast construction of ceRNA networks. SPONGE uses 'multiple sensitivity correlation', a newly defined measure for which we can estimate a distribution under a null hypothesis. SPONGE can accurately quantify the contribution of multiple miRNAs to a ceRNA interaction with a probabilistic model that addresses previously neglected confounding factors and allows fast P-value calculation, thus outperforming existing approaches. We applied SPONGE to paired miRNA and gene expression data from The Cancer Genome Atlas for studying global effects of miRNA-mediated cross-talk. Our results highlight already established and novel protein-coding and non-coding ceRNAs which could serve as biomarkers in cancer. AVAILABILITY AND IMPLEMENTATION: SPONGE is available as an R/Bioconductor package (doi: 10.18129/B9.bioc.SPONGE). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide enhancer annotations differ significantly in genomic distribution, evolution, and function.

BACKGROUND: Non-coding gene regulatory enhancers are essential to transcription in mammalian cells. As a result, a large variety of experimental and computational strategies have been developed to identify cis-regulatory enhancer sequences. Given the differences in the biological signals assayed, some variation in the enhancers identified by different methods is expected; however, the concordance of enhancers identified by different methods has not been comprehensively evaluated. This is critically needed, since in practice, most studies consider enhancers identified by only a single method. Here, we compare enhancer sets from eleven representative strategies in four biological contexts. RESULTS: All sets we evaluated overlap significantly more than expected by chance; however, there is significant dissimilarity in their genomic, evolutionary, and functional characteristics, both at the element and base-pair level, within each context. The disagreement is sufficient to influence interpretation of candidate SNPs from GWAS studies, and to lead to disparate conclusions about enhancer and disease mechanisms. Most regions identified as enhancers are supported by only one method, and we find limited evidence that regions identified by multiple methods are better candidates than those identified by a single method. As a result, we cannot recommend the use of any single enhancer identification strategy in all settings. CONCLUSIONS: Our results highlight the inherent complexity of enhancer biology and identify an important challenge to mapping the genetic architecture of complex disease. Greater appreciation of how the diverse enhancer identification strategies in use today relate to the dynamic activity of gene regulatory regions is needed to enable robust and reproducible results.

Von Hippel-Lindau Acts as a Metabolic Switch Controlling Nephron Progenitor Differentiation.

BACKGROUND: Nephron progenitors, the cell population that give rise to the functional unit of the kidney, are metabolically active and self-renew under glycolytic conditions. A switch from glycolysis to mitochondrial respiration drives these cells toward differentiation, but the mechanisms that control this switch are poorly defined. Studies have demonstrated that kidney formation is highly dependent on oxygen concentration, which is largely regulated by von Hippel-Lindau (VHL; a protein component of a ubiquitin ligase complex) and hypoxia-inducible factors (a family of transcription factors activated by hypoxia). METHODS: To explore VHL as a regulator defining nephron progenitor self-renewal versus differentiation, we bred Six2-TGCtg mice with VHLlox/lox mice to generate mice with a conditional deletion of VHL from Six2+ nephron progenitors. We used histologic, immunofluorescence, RNA sequencing, and metabolic assays to characterize kidneys from these mice and controls during development and up to postnatal day 21. RESULTS: By embryonic day 15.5, kidneys of nephron progenitor cell-specific VHL knockout mice begin to exhibit reduced maturation of nephron progenitors. Compared with controls, VHL knockout kidneys are smaller and developmentally delayed by postnatal day 1, and have about half the number of glomeruli at postnatal day 21. VHL knockout nephron progenitors also exhibit persistent Six2 and Wt1 expression, as well as decreased mitochondrial respiration and prolonged reliance on glycolysis. CONCLUSIONS: Our findings identify a novel role for VHL in mediating nephron progenitor differentiation through metabolic regulation, and suggest that VHL is required for normal kidney development.